Enhanced Weiss oscillations in graphene 
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The magneto-conductivity of a single graphene layer where the electrons are described by the Dirac 
Hamiltonian weakly modulated by a periodic potential is calculated. It is shown that Weiss oscilla- 
tions periodic in the inverse magnetic field appear, that are more pronounced and less damped with 
the increment of temperature as compared with the same oscillations in a typical two-dimensional 
electron system with a standard parabolic energy spectrum. 
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I. INTRODUCTION 

Recently the successful preparation of monolayer 
graphene films 1,2 : 3 has generated a lot of activity in the 
physics of two-dimensional (2D) Dirac fermions. The 
massless energy spectrum and the specific density of 
states of electrons and holes described by the Dirac 
Hamiltonian enabled to study experimentally the chiral 
tunnelling and Klein paradox in graphene^, and led to 
the discovery of the unconventional "half-integer quan- 
tum Hall effect"— The presence of holes in graphene 
with 2D Dirac-like spectrum was confirmed by measure- 
ments of de Haas-van Alphen and Shubnikov-de Haas 
(SdH) oscillations Z These magnetic oscillations appear 
due to the interplay of the quantum Landau levels with 
the Fermi energy in the metal, and serve as a powerful 
technique to investigate the Fermi surface and the spec- 
trum of electron excitations. 

Another technique which was successfully used to ob- 
tain information on the electron spectrum of 2D systems 
is based on the interaction of electrons with artificially 
created periodic potentials with periods in the submi- 
crometer range. Such electrical modulation of the 2D 
system was created by two interfering laser beams^, or 
by depositing an array of parallel metallic strips on the 
layer surface^, and led to the discovery of Weiss oscilla- 
tions in the magneto-resistance. These oscillations are a 
consequence of the commensurability of the electron cy- 
clotron orbit radius at the Fermi energy and the period 
of the above electrical modulation. They were found to 
be periodic in the inverse magnetic field like the SdH 
oscillations, but have a different period versus electron 
density dependence. The period for Weiss oscillations 
varies with electron density (n e ) as yjn~ e , whereas that of 
the SdH ones as n e . Theoretical calculations of these os- 
cillations were presented in Refs^ ' 10 : 11 , and it was shown 
that Weiss oscillations in the magneto-resistance for mo- 
tion perpendicular to the oscillating potential can be un- 
derstood as being a semiclassical effecti 12 : 13 

The Klein paradox in graphene where Dirac electrons 
can penetrate through potential barriers with a rather 
high probability shows that electric control of Dirac elec- 



trons cannot be realized. Therefore, it is interesting to in- 
vestigate the sensitivity of Dirac electrons on the electri- 
cal modulation of the layer. Thus, we subjected the sys- 
tem to a periodic potential that introduces a new length 
scale and a new energy scale into the problem. We found 
that such a periodic potential also leads to Weiss oscilla- 
tions in graphene which are even more pronounced than 
in typically 2D electron gases with a parabolic electron 
spectrum. 

The paper is organized as follows. In the next Sec. all 
necessary expressions for the magneto-conductivity cal- 
culation are given, and in Sec. Ill the obtained results for 
graphene layer are compared with results for the standard 
2D electron system. In Sec. IV the asymptotic expression 
valid in the quasi-classical region is obtained. In Sec. V 
some simple classical explanation of obtained results are 
presented, and in the last Sec. the short conclusions are 
given. 



II. ELECTRICAL MAGNETOTRANSPORT 

We consider the graphene layer within the single elec- 
tron approximation where the low energy excitations 
are described by the two-dimensional (2D) Dirac-like 
hamiltonian* 



Hq = v F cr (^-ihV + -A 



(1) 



Here er = {(T x ,a y } are the Pauli matrices, and the vector 
potential A — {0, Bx} describing the magnetic field B = 
{0, 0, B} perpendicular to the graphene layer is chosen in 
the Landau gauge. The parameter vf characterizes the 
electron velocity which is usually about 300 times smaller 
than the velocity of light. The total Hamiltonian consists 
of two parts 



H = H + U(x), 
where the additional potential 

U(x) = V cos(27rx/a ) 



(2) 



(3) 



describes the static electrical modulation of our 2D sys- 
tem in the x direction. 
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The conductivity tensor a^ v (uj) within the one-electron 
approximation was evaluated in Rcf. 14 . We shall restrict 
ourselves to the diffusion contribution (i. e. which is the 
dominant contribution) which stems from the diagonal 
part of the density operator. Following Rcf. 15 it can be 
written as 



where 



(3e 2 

L X Ly 



£/(#c){i 

C 



f(E c )}r(E c )(vtf. (4) 



Symbols L Xl L y characterize the dimensions of the layer, 
and (3 is the inverse temperature. This diagonal part is 
caused by the influence of the electrical modulation on 
the electron drift in crossed electric and magnetic fields 
and is in order of magnitude larger than the other diag- 
onal component a xx which appears due to scattering on 
imperfections. Symbol £ denotes the quantum numbers 
of the electron eigenstates, and the velocities v% can be 
calculated as derivatives of the energy eigenvalue over the 
corresponding electron momenta. 

The resistivity tensor is given in terms of the in- 
verse conductivity tensor, namely, p xx = a yy / \a xx a yy — 
a xy a yx ) what in our case of small conductivity reduces to 
p xx « &yy/& xy , and its p xx component is actually propor- 
tional to the conductivity in the considered perpendicular 
direction. 

The energy eigenvalues are defined through the solu- 
tion of the stationary Schrodinger equation 



{H-E}H!(r) = 



(5) 



with total Hamiltonian We will follow RefJ^ and 
assume that the electrical modulation is small enough 
and restrict to a lowest perturbation expansion in Vo- For 
this purpose we have to solve the zero order Schrodinger 
equation with the Hamiltonian Hq. Due to the system 
homogeneity along y axis we substitute the eigenfunction 
as 



1 



Jk y y 



a{x) 
b(x) 



(6) 



and transform the zero order eigenvalue problem into the 
following two ordinary differential equation set for the 
wave function components: 







d 



ihv F i lr -i ir +r,)b-Ea = 
ox ay I 



(7a) 
(7b) 



where / = y/ch/eB is the magnetic length. Solution of 
these equations can be easily obtained making use of the 
analogy with the harmonic oscillator eigenvalue problem. 
It reads 



25? = V -^^ 



e -x 2 /2 

$ n (x) = H n (x), 



(9) 



is expressed in the normalized Hcrmitian polynomials, 
and x = l 2 k y indicates the localization of the electron 
in the x direction. 

For the first order correction one has to add to the en- 
ergy eigenvalue the diagonal matrix element of potential 
([3j> calculated with the above eigenfunctions: 



AE ntky = dx dyV+ k (r)U(x)y ntkv {r) 

J-oo Jo (10) 

= Yl C os(Kx )e- u / 2 {L n (u) + L n _i(u)} , 

where K = 27r/ao, u = K 2 l 2 /2, and L n (u) is a Laguerre 
polynomial. This energy correction makes the degenerate 
Landau levels ^-dependent, expands them into bands, 
what finally leads to non zero velocities: 



v y = = !-tA£„ 



h dk 



T^ue~ u/2 [L n (u) + L„_i(u)] sm(iTa;o). 
riK 



(11) 



Now substituting the velocities in Eq. |QJ and specify- 
ing the summation over quantum numbers as 



T r L m /l 2 oo 



C 



(12) 



n=0 



we obtain the final expression for the diffusion contribu- 
tion for the dc conductivity under the consideration 



where 

and the function 
4> 



A = 8tt 2 



e 2 V 2 rP 
h h ' 



(13) 



(14) 



kT [f /''"L [L n (u) + L n ^(u)] 2 (15) 

will be considered as a dimensionless conductivity. In or- 
der to simplify this expression we introduced the follow- 
ing exponential function f(E) — exp{(3(E — Ep)} where 
Ef is the Fermi energy. 



III. COMPARISON WITH THE USUAL WEISS 
OSCILLATIONS 



e ik y y 



y/2Lyl 



-i$ n -i([x + Xq\/1)\ , 
$ n ([x + xo]/t) J' [ j 



It is interesting to compare the obtained expression for 
dc conductivity with the same conductivity calculated in 
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RefJ^ for the case of a 2D electron system localized at 
the interface between two semiconductors and having a 
parabolic energy dependence on its momentum. Com- 
paring velocity expression (jlip with the equivalent one 
in ReLi 3 . (see Eq. (4)) we notice that the diffusive con- 
ductivity for the system of Dirac electrons differs from the 
one for the standard interface case by replacement of the 
Laguerre polynomial L n (u) by the average of two succes- 
sive polynomials [L n (u) + L n -x(u)] /2. And of course, 
two different expressions for the Landau level energies 
(Eq. I|8ap for Dirac electrons and 



eB 

E„ = fhj c (n + 1/2), lj c = — 

mc 



(16) 



for the parabolic electron spectrum case) has to be used. 

These differences in the expressions, however, lead to 
essentially different results for the dimensionless conduc- 
tivity as shown in Fig. [1] The results are shown as func- 
tion of the inverse magnetic field for the temperature 
T = 6K, electron density n e — 3 • 10 n cm~ 2 , and the 
period of electric modulation ao = 350 ran. The dimen- 
sionless magnetic field b = B/Bq is introduced using the 
characteristic magnetic field B = chjea\ corresponding 
to the magnetic length equal to the modulation period 
ao, which in the above case is equal to Bq = IT. 



e 




B„/B 



FIG. 1: The dimensionless conductivity versus inverse mag- 
netic field: 1 - Dirac electrons, 2- electrons with the parabolic 
energy spectrum. 



We see that in graphene (curve 1) the Weiss oscilla- 
tions are more pronounced as compared with the system 
of electrons with the standard parabolic energy spec- 
trum (curve 2). Furthermore, the Weiss oscillations in 
graphene are much more robust with respect to temper- 
ature damping in the quasi-classical region of small mag- 
netic fields. The physical reasons for these differences lay 
in different Fermi velocities of Dirac and standard elec- 
trons. In order to confirm the above statement we shall 
consider the asymptotic behavior of Weiss oscillations in 



the quasi-classical region which according to Refs^^ de- 
scribes the main features of the above oscillations and 
which also allow for explicit analytical expressions. 



IV. ASYMPTOTIC EXPRESSIONS 

The asymptotic expression for the conductivity (|15| is 
obtained following the approach of Ref.— for case of stan- 
dard electrons. That approach is applicable when many 
Landau levels are filled, and it is based on the following 
asymptotic expression for the Laguerre polynomials: 



e-^L n {u) -> -== 
Taking the continuum limit 



cos (2^/nu - tt/4) . (17) 
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and having in mind that u = 27r 2 /6, we transform 
Eq. (fT5|) into the following integral: 



I 

vph 



f(E)EdE 
o [f(E) + l] 2 ^ 



f(E)dE 
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2\/nu — — 



irvph 



Ea 



(19) 



I vpbh 



Now assuming that the temperature is low (fi~ x <C 
Ep) and replacing E = Ep + we rewrite the above 
integral as 



^ 2a ° 

q> = COS — , 

v F hbP \p J J_ oa (e s + l) 2 

r, ( 2itp ir 2np 
x cos 2 -r- - 7 + -r§s 
\ b 4 bp 



(20) 



where symbol 



p = — ^— - = kpdQ = \/2nn e ao 
vfTi 



(21) 



stands for the dimensionless Fermi momentum of the 
electron. Note in Eq. (|20|) we replaced all energies E by 
the Fermi energy Ep except that one which is included in 
the last cosine function, where the small energy correc- 
tion can influence the damping of the Weiss oscillations. 

The obtained expression for the dimensionless conduc- 
tivity can be calculated using the standard integral 



cos(air) 

cosh(f3x) ~ 2/3sinh(a7r/2/3)' 



dx 



(22) 
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and the result can be presented as 



be presented as follows: 



$ = 




where 



T 

Air z l£, \p 
+ 2A ( — ) cos 2 



A(x) 



I- A 



T 
T~r> 



2n 



P 1 



(23) 



sinh(x) 



= 2xe 



and the symbol Tjj defined as 



kpTp, 



bh 



47r 2 a 



v F 



(24) 



(25) 



(ks is the Boltzman constant) gives the temperature 
scale for damping of the Weiss oscillations. 

The validity of the asymptotic expression is seen in 
Fig. [2] where the Weiss oscillations for the Dirac elec- 
trons (curve 1) are shown together with their asymptotic 
expression (curve 2) for the same parameter values as in 
Fig. [1] in the region of strong magnetic fields where devi- 
ations of the asymptotic expression to the exact expres- 
sion are largest. We see that the coincidence of the exact 



T = 6K 

n a = 3 10" cm 

a„ = 350 nm 




B/B 



T 



4ir 2 T P 




A 



2-1 I jr ) <-<>* 2 



T 



b 8 



(26) 



where the critical temperature reads 

bph 2 



knTp = 



4n 2 i 



Having in mind that 



P_ 

m 



kpiio 



a 



F i 



(27) 



(28) 



where Vp is the velocity of the standard electron on the 
Fermi surface, the critical temperature can be presented 
as 



k B T P = 



bh 



47r 2 ao 



F ■ 



(29) 



For parameters used in Fig.[T]plot kp ~ 1.4 • 10 6 , p ~ 50, 
and cos(7r/p) ~ 1. Thus the asymptotic behavior of the 
dimensionless conductivity for Dirac electrons J23|) and 
standard electrons with parabolic energy spectrum J26|) 
differ mostly due to the very different critical tempera- 
tures. 

Comparing Eqs. (|25]) and (|29|) we see that 



Tp 



U F 
Vp' 



(30) 



or the ratio of critical temperatures for the standard and 
Dirac electrons is equal to the ratio of the corresponding 
velocities on the Fermi surface. It can be estimated as 



Tp h\J2nn e 
Tp, vpmaQ 



(31) 



and for typical cases is less than unity. For instance, in 
the case of the parameters used in Fig. |T] it is Tp/Tp, ~ 
0.24 what explains the different slope and damping of the 
Weiss oscillations. 



FIG. 2: The graphene conductivity as a function of the in- 
verse magnetic field: 1 - exact solution (|15|l , 2 - asymptotic 
expression (|23|l . 



result and its asymptotic expression is rather good ev- 
erywhere except the region of very strong magnetic field 
where SdH oscillations become superimposed on top of 
the Weiss oscillations. 

We compare now the obtained result for the conductiv- 
ity for the Dirac electron system with the similar asymp- 
totic result for the system of electrons with the standard 
parabolic energy spectrum taken from RefJ^ which can 



V. QUASI-CLASSICAL EXPLANATION 

The obtained results can be understood from a simple 
physical picture. In order to estimate the oscillation pe- 
riod we write down the momentum of the electrons on 
the Fermi surface: 

pp = muj c R c (32) 

which is identical in both cases. Then it follows that the 
radius of electron orbit in the magnetic field is 

pp hk F c pi 2 

li c = = — — = . (33) 

mojc en ao 
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The physical reason for the appearance of the Weiss os- 
cillations is the commensurability of the electron orbit 
radius with the period of the electrical modulation, con- 
sequently, the argument of the cosine function has to be 
proportional to 



Rc 
«o 



I 2 



P- 



P 
b' 



(34) 



dividing Eq. ([36)1 by Eq. ([37)) we obtain the same critical 
temperature ratio dependence on the ratio of the Fermi 
velocities ([30]) what confirms the quasi-classical nature of 
Weiss oscillations. 



VI. CONCLUSIONS 



which we can see in both expressions (|23[) and ([28)) . 

The damping factor of the oscillations can be estimated 
as follows. Due to the finite temperature there are elec- 
tron orbits with various radia. The effective damping can 
be estimated as the ratio of the above broadening of the 
orbit to the period of the modulation: 



ao oEf 



(35) 



In the case of standard electrons (when Ep = p F /2r 



it reads 



dR r 



a fip F dp F 
k B T 



ao/3hkpmu> c 



ma,Q 



(36) 



fc R r, 



in agreement with the definition of the critical tempera- 
ture for standard electron ([27)) . 

In the case of Dirac electrons (when Ep = vfPf) the 
above damping parameter can be estimated as 



Sd = 



1 dR r 



ck B T 



fiaoVF dpF aoVFeB HbvF 



k B T, 



(37) 



what coincides with the obtained earlier result for the 
critical temperature up to the same constant 4-7T 2 . Now 



In conclusion, we studied the Weiss oscillations in elec- 
trically modulated single layer of graphenc. It was shown 
that the static conductivity oscillations are periodic in 
1/B with period varying with electron density as y/n^ 
like in the 2D electron system with standard parabolic 
energy spectrum. Due to the larger Fermi velocity the 
conductivity oscillations of the Dirac electron system are 
more pronounced and less damped as compared with the 
2D system of electrons with parabolic energy spectrum in 
case of analogous parameters. The found Dirac electron 
sensitivity to the electric perturbation does not contra- 
dicts the Klein paradox, because in contrary to electron 
tunneling through barriers where both the electron and 
the hole nature of the excitation plays a role, only the 
electrons at the Fermi energy are responsible for the con- 
ductivity and for the studied Weiss oscillations. Thus, 
their behavior in the electric field is not overshadowed 
by the admixture of the hole states. 
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